Polarons with a twist 
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We consider a polaron model where molecular rotations are important. Here, the usual hopping 
between neighboring sites is affected directly by the electron-phonon interaction via a twist-dependent 
hopping amplitude. This model may be of relevance for electronic transport in complex molecules 
and polymers with torsional degrees of freedom, such as DNA, as well as in molecular electronics 
experiments where molecular twist motion is significant. We use a tight-binding representation and 
find that very different polaronic properties are already exhibited by a two-site model - these are 
due to the nonlinearity of the restoring force of the twist excitations, and of the electron-phonon 
interaction in the model. In the adiabatic regime, where electrons move in a toiu-frequency field 
of twisting-phonons, the effective splitting of the energy levels increases with coupling strength. 
The bandwidth in a long chain shows a power-law suppression with coupling, unlike the typical 
exponential dependence due to linear phonons. 
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There has been considerable interest in the study of 
electronic transport in organic materials. They provide 
technologically relevant alternatives for new electronic 
and optical materials, and offer an opportunity to gain 
deeper understanding of quantum transport processes. 
Moreover, transport in molecules also constitutes the ba- 
sis for important biological phenomena, such as photo- 
synthesis and other charge transfer eventsE Recent de- 
velopments in the fabrication of high electron mobility 
systems using organic semiconductors, have further em- 
phasized the need to understand mobile polaronic behav- 
ior in complex molecular crystals.EIln addition, a series of 
beautiful experiments in the area of molecular electron- 
ics have underscored the importance of atomic vibrations 
in the njejasj|H-ed current-voltage characteristics of these 
systems BqeIeI 

In typical inorganic materials the atomic degrees of 
freedom (DOF) associated with translation are the most 
important in electronic transport. In contrast, the twist- 
ing DOF (or "librons" ) of molecular complexes may also 
play important roles in organic and molecular solids. 
In a wide variety of organic linear-chain crystals of 
stacked molecules, the librons may in fact be essen- 
tial. For instance, it was found that libron scatter- 
ing was important in the understanding |o£,|CpBductiv- 
ity mechanisms in the TTF-TCNQ system^yBllias well 
as in phthalocyanine-based "molecular metals" .M More 
recently, the conductivity of DNA and (pnsti±JiPi|it .base 
crystals has beeaiSludied experimentallyJ13ll3oHE3 and 
theoretically.lllHlall3'Ea^El Since DNA can be viewed as 
a "molecular wire" , a one-dimensional twisted chain pi 
stacked base pairs with a somewhat flexible structure,E3 
it is likely that parallel rotation of adjacent base pairs 
will be important in electronic transfer mechanisms. No- 
tice that these modes would have a stronger impact on 
the electronic overlaps along the chain than other flex- 
ural modes. In fact, a variable range hopping theory 



where the twist angle fluctuations between base pairs are 
taken into account was used recently to explain the tem- 
perature dependence, of electrical conductivity along the 
DNA double helix. E£l Bruinsma et al. have also studied 
charge migration along DNA based on a model in which 
large structure fluctuations, including [Relative rotation 
between adjacent bases, are considered^ 

There exist a few models of electron-libron interactions 
in the literature. Many of them have considered the 
situation of small rotation ajjigles.in. a. harmonic or nearly 
harmonic approximation.BE£lEJ'Ej£3 Although-j-other 
authors have considered large angle rotations,t§E3 the 
molecular kinetic energy term was neglected, and thus 
the inertial backreaction effects in molecules were not 
taken into account. Our model considers the possibility 
of large rotation effects, which are not negligible in 
systems such as (TTF)7l5, and others with similarly 
flexible structure. We also include the full kinetics of the 
molecule in order to include backreaction effects. Our 
description is fully quantum mechanical and includes 
strong anharmonic interactions, shown to be essential 
for polaron self-trappingo From the full solution of a 
two-site model, we flnd that in the adiabatic regime, 
where electrons move in a low frequency field of twisting- 
phonons, the splitting of the energy levels (and therefore 
the effective intersite hopping) increases with phonon 
coupling strength. In the opposite, non-adiabatic 
regime, the effective hopping constant increases with 
phonon frequency, but saturates at high values. For a 
long chain model with local- twists, we find a power-law 
suppression of the bandwidth, unlike the exponential 
result for translational polarons. These results arise 
from the different coupling scheme and the nonlinear 
twist-polaron restoring force, as we now describe. 

The tvi^ist polaron model. The on-site interaction of 
the electron with phonon modes tends to localize the elec- 
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tron. On the other hand, the molecular transfer (hop- 
ping) integral in a given atomic or molecular configura- 
tion results in electron delocalization. In an extended 
system, the hopping integral competes with the phonon 
interaction, and depending on the relative strengths of 
the different terms, one obtains a mobile polaron (wi 
heavier effective mass), or a self- localized polaron. I 
Our model allows for such physical behavior as well. 
Most importantly, however, the phonon enters here mod- 
ifying also the hopping integral and is then substan- 
tially different to typical models for translation-phonon 
polarons.EZi 

For simplicity, let us first consider a two-site model. 
The dynamics of the system is described by the Hamil- 
tonian 



Hip = E'V', we use 
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H = -(c+ci 

-{to + ti cos6){clc2 + cjci) 

21 de^ 



(1) 



K cos 9 , 



where the first term represents the possible on-site energy 
modulation due to the molecular rotations with angular 
DOF 9, and the cf operators describe the two electronic 
wave functions ( "sites" ) . The second term is the rotation- 
dependent hopping term. Notice that the molecular ro- 
tation would either enhance or suppress electron hop- 
ping in this model, depending on whether ti and Iq have 
the same or opposite sign. A microscopic calculation of 
stacked molecules, such as the base pairs in DNA, or the 
phthalocyanine stacks, allows one to estimate the effept 
of molecular twists on the electron hooping or overlap.E3 
Recent results for guanosine crystals,E3 for example, show 
clearly that eclipsing of molecules in successive "layers" is 
essential for substantial electron hopping between them. 
This twist dependence is introduced in our model via 
the hopping amplitude ii, and then plays the role of a 
coupling constant. The last two terms in (1) are the rel- 
ative kinetic energy for the molecules and the interaction 
potential characterized by the libron restoring 'force con- 
stant' K. Here, 9 = 9i—92, where 9i is the rotation angle 
of each molecule, and / is the reduced moment of iner- 
tia for the relative rotation of the two adjacent molecules 
(see inset Fig. 1). An overall rotation of the pair is gap- 
less and neglected. Notice also that the first term in (1) 
leads to a constant shift of K, which is then implicit in 
the following discussion. Notice that this model consid- 
ers nonlinear coupling to the electrons via ii,Lfl as well 
as a non-linear restoring force for the relative rotation of 
this two-molecule stack. On the other hand, the model 
does not consider the possible mixing of rotation DOF 
with translation displacements, although this coupling is 
for the most part smalLEj We believe that our model 
contains the main physics of polaronic excitations corre- 
sponding to the rotation DOF. 

In the following, we use the energy unit for con- 

venience, so that all energy quantities are dimensionless 
from now on. In order to solve the eigenvalue problem 



6(9)4 \o) + m4\o) 



(2) 



which yields, from {H — el)tp — 0, a set of secular equa- 
tions for S and (3. Defining A = S + (3, and B = 5 — f3, 
we have 



{-^ +e + to+pcos9)A = 
i-^+e-to + qcos9)B = 0, 



(3) 



with q = K — ti, and p + ti . The solutions are peri- 
odic Mathieu functionsJ2j with properties which depend 
on the values of the coefficients p and q. 

As we will see below, the rotation DOF in the model 
produce modifications of the electronic hopping. In 
order to reveal the physics more clearly, we analyze the 
two-site problem in more detail in the adiabatic and 
non-adiabatic regimes. 

Non-adiabatic regime. In the non-adiabatic regime 
K 3> tQ,ti, the electrons move in a field of high frequency 
twist-phonons. We first neglect fully the hopping terms. 
We find two degenerate states 



^1 = 







V'2 = 







(4) 



where (f>o is the normalized ground state (with eigenvalue 
Sq), solution of the equation 



d^4 
'dp 



[e + K cos 6] 



0. 



(5) 



To calculate the first order correction, we take the wave 
function = aipi + /3V'2, so that the secular equation 
becomes 



= 0, 



(6) 



where the renormalized twist-dependent hopping con- 
stant is 



t[ =tij' 



cos 9 (j)nd9 . 



(7) 



The eigenvalues are simply £i_2 = £o ± \to + t[\, so that 
the splitting of the two originally degenerate eigenvalues 
is Ae = 2\to + t[\. It is easy to see from Eq. (0) that 
t'l < h- Therefore, Ae < T, where T = 2to + 2ti is the 
hopping constant "without rotation" , i.e. = 0. 

In the regime of large K {K ^ 1), the local twist 
potential is effectively harmonic and one can write the 
solution as a superposition of local harmonic oscillators 
about 2m7r, 



a 



1/2 



H ^1/4' 



-Q^(e-2mir)^/2 



(8) 
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where a"^ = K/2, and the eigenvalue is Eq 
Correspondingly, 



-K 



(9) 



and Ae = 2to + 2ti{l-l/4a'^) < T, with Ae approaching 
T in the large K limit. In this case, the molecule is weakly 
oscillating about 9 — (mod 27r) , with a relatively minor 
impact on the electron hopping. Notice that the resulting 
"bandwidth narrowing" is similac-in form to the result 
of a two-site vibrational polaron.cZl There, the effective 
bandwidth in the non-adiabatic regime is J = Je~'^ 
where J is the original bandwidth and 7 is the electron- 
phonon coupling constant. Notice however the different 
role that K plays here, and that the bandwidth is only 
weakly suppressed in this case. 

When K is small (1 ^ toiti), the ground state 

function can be written asEB 



1 



v/27r(l + i^V2) 



(1 + /rcos^ 



(10) 



and the eigenvalue is eo = ^^I"^^- Wc find that t[ — 

h-, — ^TTTT - Kti, and Ae = 2to + 2Kti < T. Notice 
l + K^/2 

that in this limit of soft restoring force, the molecules 
are nearly uniformly rotating at all angles, and strongly 
suppress the electron intersite hopping. 

We see that with increasing K, the renormalized 
twist-dependent hopping constant t[ increases and 
saturates to ti, as shown schematically in Fig. 1. We 
should emphasize that one can view the resulting Ae 
as akin to a bandwidth in the long chain limit. These 
results then indicate that the twisting DOF strongly 
prevent the electron from propagating along the chain. 
The "stiffer" the rotation mode or libron, the easier for 
the electrons to move from site to site. 

Adiabatic regime. In this limit, K <C to,ti, the elec- 
tron moves 'rapidly' in a field of low frequency 'slow' 
phonons. We assume a wave function of the form 



ai9) 

m 



3(6) 



(11) 



and first consider the electron moving in a frozen back- 
ground. 



0. 



-Eo{d) -{to + hcos9)\ f si{9) 
-{t„+hcos0) -En{9) )\S2{9) 

(12) 

The corresponding eigenvalues — ±\to + ti cos 6*1 play 
the role of a potential energy term in the equations for 
3(9), 



89^ 



Kcos9-E^{9)+e}S{9) = 0. 



(13) 



These equations can also be solved in terms of Mathieu 
functions. It is easy to find that when K ^ to,ti, and 
ti ^ 1, the eigenvalues of the problem are given by 



El 



£2 



-op ''O 



-op 



to + K- 




(14) 



where Eop — —ti + \ft\j2 is the polaronic energy shift 
due to the electron-libron off-site interaction. In this 
case, Ae = 2to + 2K — K j \/2t\. Unlike the non-adiabatic 
case, Ae saturates to 2io + 2K for increasing ii, and 
it therefore shows a strong bandwidth narrowing. 
Notice, however, that the two-site adiabatic problem 
for the small Holstein polaron yields an-|even stronger 
suppression of tunneling. In that case,EII one finds an 
exponential narrowing with coupling, ~ e ''' similar 
to the non-adiabatic case mentioned above, but with 
large 7 values. 



Long chain rotational polaron. The long-chain gen- 
eralization of the two-site model is an interesting phys- 
ical problem with additional matiiematical complexity 
which will be reported elsewhere. lj Here, we focus on 
a simpler local-twist model involving local angular coor- 
dinates. The local-twist hamiltonian is given by 



H = 



21W 



- K cos 9i - sicla cos 9^ - to{cfci+i + c. 



(15) 



which includes anharmonic restoring force and local non- 
linear electron-phonon coupling via ei (the more general 
ti-dependence is also discussed in Ref. |o|). 

In the non-adiabatic regime, K ^ 1, one can gen- 
eralize the two-site approach and identify two limiting 



cases. For strong restoring force, K ^ Si, it is possible 
to show that the polaron energy shift is — ei, while the 
bandwidth is slightly reduced to A = 2to(l — Ei/K). 
In the case of softer twist frequency, ei ^ iiT ^ 1, 
the bandwidth is given by 2to{K/eiY/^ . This power 
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law suppression of hopping is much weaker than the. 
usual Holstein 'translational' polaron in long chains.EZI 
The different dependence arises mostly from the local 
anharmonic coupling in this case, and would result in 
the polaron state having a higher mobility along the 
chain, despite the strong coupling constant. 

Conclusion. We have introduced the idea of librons in- 
teracting with electrons in a complex molecular system. 
The strong coupling of rotational DOF yields a unique 
class of twisting polarons. Inclusion of anharmonic 
coupling to the electron via hopping and a nonlinear 
restoring potential, were shown to result in rather dif- 
ferent polaronic behavior. We find electronic bandwidth 
suppression thanks to the excitation of twisting phonons, 
just as in the case of the usual vibrational DOF, but with 
a comparatively weaker effect in our case. This difference 
is due to both of the nonlinear characteristics men- 
tioned, and provides a possible route for differentiation 
in experiments. We anticipate that the twist polaronic 
features could be tested by performing experiments in 
a variety of different molecules, where the restoring 
force is changed systematically, for example. We expect 
that our model will provide some insights on the low- 



temperature transport properties of interesting systems, 
such as recent molecular electtonics experiments, and 
helical polymer "soft" crystals .Eil Notice however that a 
detailed comparison with experiments will require one 
to consider the coupling of these polarons to charge or 
current leads. Such effects [in vibrational polarons have 
only recently been studied,E3 and one would anticipate 
qualitative differences, which will be discussed elsewhere. 
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FG02-91ER45334, NSF NIRT grant no. 0103034, and 
from the Condensed Matter and Surface Sciences Pro- 
gram at Ohio University. Helpful discussions with the 
NIRT group at OU are also appreciated. 



FIG. 1: Effective tunneling between sites, ti, as a function of 
twist mode restoring force, K. t^ is given by Eq. (9) for large 
K values, and is linear for small K, all in the non-adiabatic 
regime. Inset: diagram of parallel twist deformation of stack 
of molecular units with rotational degree of freedom along the 
stack. 
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